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Express Mail Label No. EV303487174US 



METHOD AND SYSTEM FOR 
DETERMINING DISTORTION IN A CIRCUIT IMAGE 

BACKGROUND 

[0001] Integrated circuits (ICs) are ubiquitous in today's society. ICs have 
microscopic dimensions, which may make ICs difficult to test and characterize. 
Recent techniques for measuring IC parameters include Picosecond Imaging 
Circuit Analysis (PICA). In general, PICA techniques provide a measure of the 
relative timing and source position of infrared light or photon emissions that occur 
at various locations across an IC at various times. The PICA technique relies on 
current-carrying transistors emitting infrared light as the transistors switch states. 
Because current flow is indicative of a switching transition, and because infrared 
light emission is indicative of current flow, photon detection may be used to 
measure the timing of switching events in digital circuits. 

[0002] However, accurately extracting spatial-related and timing-related 
parameters from PICA measurements is not a straightforward procedure. Photon 
emission is a statistical process, and the probability of detected photon emission 
for current very large scale integration (VLSI) technology is on the order of 10" 9 
infrared photons per switching event per micron of transistor channel width. To 
acquire a statistically useful measurement, many switching events must be 
generated, and the timing of the separately detected photons must be correlated 
with a repeatable trigger. The relative timing of the photons can be plotted as a 
histogram, and the peak of the histogram may be used to determine the timing of 
the transistor's switching transition relative to the trigger. Unfortunately, the 
histograms are "noisy" as a result of the granularity of the events (i.e., either a 
photon is detected, or not detected) and the rarity of events. The noise makes the 
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measurements difficult to interpret. The histograms gradually approach a smooth 
curve if the measurements are taken over a sufficiently long period of time, but 
this is feasible for small width devices or devices that switch infrequently relative 
to the clock frequency. Furthermore, as device sizes continue to decrease, 
histograms representing spatial photon emission will tend to become less precise. 

SUMMARY 

[0003] An objective method of extracting spatial and timing related circuit 
parameters from picosecond-scale photon timing measurements is desirable. In 
accordance with at least some embodiments, a method of determining distortion 
in a circuit image comprises: defining potential photon emission areas in the 
circuit image using a layout database, determining ideal photon emissions over 
the potential photon emission areas, measuring photon emissions for the 
potential photon emission areas, comparing the ideal photon emission with the 
measured photon emissions, and producing a mathematical model that predicts 
the amount of spatial distortion over the potential photon emission area. 

BRIEF DESCRIPTION OF THE DRAWINGS 
[0004] A better understanding of the present invention can be obtained when 
the following detailed description is considered in conjunction with the following 
drawings, in which: 

[0005] Figure 1 shows an exemplary picosecond-scale photon detection 
system; 

[0006] Figure 2 shows an exemplary functional block diagram of an 
apparatus for analyzing photon detection data; 

[0007] Figure 3 shows a PICA image of four inverters including spatial 
distortion; 

[0008] Figure 4A shows an algorithm for developing a model of spatial 
distortion; 

[0009] Figure 4B shows several Laplace based distributions with varying 
scale factors; and 

[0010] Figure 5 shows a flow diagram for determining distortion in a 
circuit image. 
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DETAILED DESCRIPTION 
[0011] Figure 1 shows an illustrative functional block diagram of an apparatus 
for making picosecond-scale photon timing measurements. The apparatus 
includes a test pattern generator 302 that generates a sequence of input signal 
vectors, or "test vectors," that are designed to make a device under test 304 (e.g., 
an IC) operate in a desired manner. The desired operation depends on various 
circumstances. For example, when a fault has been detected in the device 304 
under certain operating conditions, the test pattern generator is programmed to 
duplicate those operating conditions. When measuring propagation delay 
between circuit components, the test vectors may be designed to repeatedly 
cause transitions in those components. 

[0012] The test pattern generator 302 provides the test vectors to the device 
under test 304, which is supported in a test stand 306. The test stand 306 
accepts input signals from the test pattern generator 302 and provides the input 
signals to the device under test 304 along with power and any other needed 
signals (e.g., clock signals, bias voltages, etc.) The test stand 306 also may 
accept output signals from device 304 and provide the output signals to another 
functional block (not shown). 

[0013] Device under test 304 can be a fully packaged IC chip that has been 
modified to make the back side of the chip's silicon substrate viewable. In addition 
the substrate may be thinned to increase the transparency of the substrate to 
infrared light. 

[0014] As the device 304 operates, currents flow through various components 
within the device, and in particular, through metal-oxide-semiconductor field effect 
transistors (MOSFETs) undergoing switching transitions. The current-carrying 
components emit photons of infrared light that pass through the substrate and 
radiate from the back side. A microscope lens 308 captures and redirects the 
photons into an imaging photomultiplier 310. 

[0015] Photomultiplier 310 performs at least two functions: amplification and 
measurement. A single photon impacting on the photomultiplier 310 knocks an 
electron free from an image grid contained in the photomultiplier. The 
photomultiplier 310 amplifies the effect of the freed electron to create a cascade 
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of electrons that form a detectable current pulse. The photomultiplier 310 also 
measures the position of the incident photon in an image grid by determining 
where the current pulse was detected. The image grid may be divided into a 
suitable number (e.g., 512 or more) of rows and columns. Photomultiplier 310 
provides an X signal indicating the row, and a Y signal indicating the column in 
which the photon is detected. In addition, the transient current pulse is 
conditioned to form a T" signal that indicates the time when the photon is 
detected. The T signal provides picosecond-scale resolution. The X, Y, and T 
signals may or may not be in digital form. 

[0016] The Data Acquirer and Analyzer (DAA) 312, receives the X, Y, and T 
signals from the imaging photomultiplier 310 and receives a trigger signal from 
the test pattern generator 302. Alternatively, the trigger signal may be provided 
from the device under test 304 or the test stand 306. The trigger signal comprises 
a reference signal, and typically is correlated with the switching event(s) being 
measured on device 304. The DAA 312 measures the difference between the T 
signal relative to the trigger signal to determine photon detection time. If the 
signals are in a digital format, the DAA 312 calculates a time difference between 
the T signal and the trigger signal. If the signals are not digital, the DAA 312 
determines the photon detection time by measuring an interval between a 
transition of the trigger signal and a corresponding transition of the T signal. 
[0017] If the X, Y and T values are not in digital form, the DAA 312 may convert 
the X, Y, and T values to digital form. The DAA 312 then stores the (X,Y,T) 
triplets for later analysis. In an alternative embodiment, the DAA 312 constructs a 
histogram and then, if desired, discards the triplets. The histogram may comprise 
a three dimensional array of possible (X,Y,T) triplets, in which a count is 
associated with each of the values. When a particular (X,Y,T) triplet is received, 
the count associated with that triplet is incremented. The histogram may be 
stored in the form of a three dimensional array. Alternatively, the histogram may 
be represented using a more compact data structure — e.g., parallelpipeds as 
described below. 

[0018] The data collection performed by the DAA 312 may be divided into 
sessions. At the beginning of a session, a new data file is created. Creating a new 
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data file may include initializing the histogram with zeros. Data collection then 
commences and continues until the session is temporarily suspended or is 
ended. At the end of a session, the data collection is halted and the data file or 
histogram is exported or stored in long term storage. The data can then be 
analyzed in a variety of ways by DAA 312 or by an external system. 
[0019] Figure 2 shows an illustrative functional block diagram of a system 400 
for analyzing data obtained by the apparatus of Figure 1 . The system 400 may be 
contained within DAA 312, or alternatively, may be an external system coupled to 
DAA 312. In yet another embodiment, the system 400 is an entirely separate 
system to which the data from DAA 312 is transported. 

[0020] System 400 includes a data acquisition module 402 that receives the 
photon measurement data. Module 402 may simply be a port through which the 
data enters the system, or module 402 may be a more complex module for 
conditioning and acquiring data signals. The data is stored in a storage module 
404, which maintains a copy of the data for access by the processing module 
406. Storage module may take the form of dynamic random access memory 
(DRAM), static random access memory (SRAM), Flash ROM (read only memory), 
optical storage media, magnetic storage media, or any combination thereof. 
Storage module 404 (or a portion thereof) may be removable so as to be portable 
between systems. Storage module 404 may further maintain a copy of software 
for analyzing the data as described further below. 

[0021] Processing module 406 accepts instructions from a user via input 
module 408, and responsively operates on the data from storage module 404 to 
extract circuit timing parameters. The data or parameters may be displayed in 
some form to the user by display module 410. Storage module 404 may include 
information regarding the device under test 304, such as a layout database 411 
indicating the location of circuit elements on the device under test 304. 
[0022] The photon detection data (or a subset thereof) may be displayed and 
analyzed in a variety of forms. For example, the user may be interested in 
evaluating a transistor located at a specific X-Y coordinate during different 
intervals of time. Figure 3 depicts an extracted PICA image 500 of four serially 
connected inverters. The X and Y axes for the PICA image 500 are also shown in 
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Figure 3. The PICA image 500 indicates photon emission of the inverters over 
time. The photons generated by each of the inverters often undergo spatial 
distortion in traveling from the device 304 through the lens 308 and into the photo 
multiplier 310 so that the actual X-Y location of the photon emission is not 
accurately portrayed in the PICA image 500 that is extracted. 
[0023] Figure 3 depicts this spatial distortion in graph 502 with photon intensity 
or count from the extracted PICA image 500 on the ordinate axis and the Y pixel 
location from PICA image 500 on the abscissa axis. For ease of discussion, the X 
pixel location of each photon counted in the PICA image 500 is held constant at 
the X axis while the photons are counted at various Y pixel locations over time. 
Ideal photon emission 504 is indicated by a sharp rectangular pulse, where 
photons are emitted over an interval having a width of 2A. (The interval 2A is 
centered at Y pixel location of 1 ,675 in Figure 3, but can be scaled to 0 so that the 
interval may defined as starting at -A and ending at +A.) The ideal interval from -A 
to A corresponds to an area of a saturated transistor where uniform electric field 
strength is present across the transistor's width, resulting in idealized photon 
emissions along the area of uniform electric field strength. Measured photon 
emission 506 is illustrated in Figure 3 as a bell shaped distribution. The difference 
between the ideal photon emission 504 and the measured photon emission 506 
is referred to as "spatial smear," and may be quantified using various statistical 
methods. By quantifying the smear or distortion in the extracted data, timing 
information may be more easily extracted despite distortion. 
[0024] Based on empirical measurements akin to measurements shown in 
Figure 3, the smear in a photon's detection point can be approximated using a 
Laplace distribution with a probability density function (PDF) that has an 
exponential-power distribution. Given this PDF, the cumulative distribution 
function (CDF) may be determined by distorting (e.g., mathematically convolving), 
the ideal photon emission 504 over the interval from -A to A by the PDF. The 
coordinates for the interval from -A to A come from the layout database 41 1 . For 
example, the layout database 411 may define a transistor comprising a gate 
region where the interval of -A to A corresponds the coordinates of the gate 
region in the layout database 411. In some embodiments, points along the ideal 
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photon emission 504 may be distorted by mathematically convolving such points 
with the Laplace distribution along the interval -A to A (from the layout database 
411.) The result of such a convolution produces a modified distribution 508 that 
approximates the measured photon emission 506. 

[0025] Figure 4A illustrates an exemplary closed form algorithm 600 for 
computing the probability that a photon emitted in the interval -A to A is detected 
in the chosen time interval of (a,b), where the probability is given by the 
expression F[b,A]-F[a,A]. The expression F[b,A]-F[a,A] represents an exemplary 
closed form expression for the modified distribution 508. Closed form expressions 
are desirable because closed form expressions often require fewer computation 
cycles to achieve the desired result. Portions of the device 304 contain large 
isolated circuit elements with known photon emissions that may be used to 
calibrate the desired amount of scaling between the modified distribution 508 and 
the measured photon emission 506. For example, maximum likelihood estimation 
(MLE) techniques may be used to correlate layout database coordinates to the 
PICA image of a large isolated circuit element with known photon emission in 
order to calibrate the desired amount of scaling. The desired scaling factor, which 
is proportional to the amount of spatial smearing, may be a function of several 
fundamental factors such as backside wafer preparation, detector lens geometry, 
and the detector's spatial electrical characteristic. 

[0026] Figure 4B illustrates various Laplace based distributions generated by 
algorithms 600, where each distribution contains a different amount of scaling. 
Thus, by using information from the layout database 41 1 to find a statistical model 
that fits the distortion between the ideal photon emission 504 and the measured 
photon emission 506, a model for spatial smear may be derived. Although 
Laplace distributions are described, other embodiments may include different 
data sets where alternative PDFs are used to produce a model that fits the 
observed data. 

[0027] As an aid in making statistical inferences from PICA images, areas of 
photon emissions may be divided into "emitter" regions that are expressed as a 
set of pixel coordinates. The layout database 411 defines the possible areas of 
photon emission that may be divided into emitter regions — e.g., the channel 
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region of CMOS transistors. In general, photon emission is a rare event whose 
probability is dependent on the electrical state (i.e., electric field intensity and 
density of charge carriers of the local region). Thus, for each emitter region E, 
photon emission may be described using a Non-Homogenous Poisson Process 
(NHPP) N E (X,Y,T), where the intensity function of the emitter A E (X,Y,T) is equal to 
A E (T) if the photon intensity function is constant over its pixel coordinates but not 
constant over time. Moreover, the emitter regions in the PICA image can be 
assumed to be mutually independent processes that interact solely through their 
intensity functions with respect to time A E (T). In this manner, when modeling a 
continuous pixel region whose intensity varies, the pixel region may be divided 
into separate emitters E 0 through E N with varying intensities A 0 (T) though A N (T), 
and therefore the accuracy of the general photon intensity function A(X,Y,T) may 
be modeled to any desired level of accuracy by varying the number of contiguous 
emitters. 

[0028] Referring again to Figure 3, the measured photon emission 506 near the 
fringes of the abscissa axis (e.g., around Y pixel coordinates 500 and 2700) is 
non-zero. This non-zero measurement is related to background photon processes 
and is sometimes referred to as "dark current." Akin to modeling other photon 
emission, dark current may be modeled using an emitter with a pixel set equal to 
the entire PICA image coordinates. Furthermore, a collection of emitters may be 
used to model other background phenomena that cause gradients in photon 
emission across the PICA image. 

[0029] In addition to spatial smearing, PICA images undergo timing based 
distortion as disclosed in commonly owned patent application entitled "Method of 
Extracting Circuit Timing Parameters Using Picosecond-Scale Photon Timing 
Measurements," application Serial No. 10/269,609, filed on October 11, 2002, 
which is hereby incorporated by reference. The time distortion model developed 
in previously filed application number 10/269,609 may be generalized to be that 
of an arbitrary mixture of Gaussian distributions scaled to have a mean of zero. 
Using symbolic representation, if (j>((s-|x)/a)/a represents the PDF of a Gaussian 
time distortion model, then Equation 1 represents the PDF of a mean-zero 

mixture of such models when Y* b. =1 and Y* b.u. = 0. 



114935.04/2162.18400 



200312064 



9 



VW = T M b jAs-»j)/°j)l<rj 0) 
[0030] This general time distortion model is used to model the composite time 
jitter caused by the photon detector and its timing electronics, even when their 
time distortion distributions are highly non-Gaussian. The generalization of 
Equation 1, which uses a mixture of mean-zero approximations, advantageously 
allows closed form expressions (akin to the NHPP likelihood) that have 
computational advantages. 

[0031] The NHPP model previously developed in application number 
10/269,609 for timing distortion, along with the generalization of Equation 1 may 
be combined with the spatial distortion model described herein to produce a 
photon intensity model that accounts for both time and space distortion as given 
by Equation 2. 

^(X,Y,T)= J^ E (t-s).^(s)ds.f(x J E x ).f(y,E y ) (2) 

— oo 

The expression A E (t-s)i|j(s) represents the time distortion model (developed in 
previously filed application number 10/269,609) for an emitter region E. The 
expression f(x, E x ) represents the PDF in the X direction in the PICA image, 
where X is provided by the PICA image and the term E x is provided from the 
layout database 41 1 . Likewise, the expression f(y, E y ) represents the PDF in the 
Y direction in the PICA image, where Y is provided by the PICA image and the 
term E y is provided from the layout database 411. The time distortion model, the 
spatial distortion model in the X direction, and the spatial distortion model in the Y 
direction are each described using NHPP techniques and are therefore 
independent of each other. Due this independence, a form of likelihood 
estimation may be used to develop a closed form composite model where the 
time distortion model is "weighted" by the X and Y spatial distortion models. 
[0032] In some embodiments, computational efficiencies may be obtained by 
realizing that over time, some X-Y coordinates do not change in photon intensity. 
For example, the data triplets (X,Y,T) may be parallelpipeds such that each of the 
recorded data vectors X, Y, and T may be different lengths because the photon 
emissions at a particular X-Y coordinate remain constant from one time interval to 
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the next and may be represented once in a vector. In this manner, fewer 
mathematical calculations may be needed since vector lengths are shorter. 
[0033] Figure 5 illustrates a flow diagram of an algorithm 700 for determining 
distortion in a circuit image such as the circuit image shown in Figure 3. In block 
702 a potential photon emission area is defined using the layout database 411. 
This block may include defining a potential photon emission area for each 
background phenomena in the circuit image. Block 704 includes modeling the 
dark current of the circuit image. Photon emissions are measured for a potential 
photon emission area in block 706. Per block 708, the measured photon 
emissions are represented with unequal vectors, which may reduce mathematical 
computations as described above. Block 710 includes determining the expected 
level of photon emissions over the potential photon emission areas defined by the 
layout database. In block 712, a composite time-space distortion model is formed 
by weighting a spatial-distortion model by a time-distortion model. In some 
embodiments, the time and space distortion models are mutually exclusive and 
may include evaluating Equation 2 shown above. In block 714, the expected 
photon emission level is compared with the measured photon emission level. Per 
block 716, the amount of spatial distortion is predicted for each potential photon 
emission area based on the results of comparing the expected and measured 
photon emission levels. This predicting may include implementing a PDF or 
mathematically convolving measured photon emission levels with a PDF. 
[0034] Implications of the embodiments disclosed are numerous. For example, 
photon emissions relate to the switching of transistors and the PICA images may 
be used to determine if a transistor is not functioning properly. Because each 
successive generation of devices are smaller and smaller, transistors become 
more congested thereby increasing the difficultly in determining a transistor's 
functionality using PICA images. However, by reducing the timing and spatial 
distortion in PICA images, sharper PICA images may be produced, thereby 
increasing PICA'S usefulness as device sizes shrink. 
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